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Abstract 

Employing the Lagrange inverting series, a solution of a mixed expo¬ 
nential/quadratic transcendental equation, that can be considered a gen¬ 
eralization of the equation defining Lambert W-function, has been found 
in terms of Bessel orthogonal polynomials, and a novel representation for 
Bessel polynomials has been found. Some other cases of transcendental 
equations formally solved by means of classical orthogonal polynomials 
have been discussed, suggesting a link between Rodrigues formulas and 
the terms of Lagrange series to be investigated. 


1 Lambert-like transcendental equation 

We consider the following transcendental equation: 

(x — a)(x — b) = le x (1) 


This transcendental equation emerges in the study of the electron in the po¬ 
tential of Upon (see in) , and it can be considered a generalization of Lambert 
W-function defining equation. Indeed, in the limit case -b= 1 =oowe obtain: 

x — a — le x 

that can be solved by means of Lambert W function: 

x = a — W{—le a ) 

Now we will produce a solution, based on Lagrange inversion series, of equa¬ 
tion [T| 

We can rewrite such equation as follows: 

le x 

x = a H-- 

x — o 
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where only trivial algebraic manipulations are employed. 
Now, remembering that the solution of the equation : 

x = a + lf(x) 

can be written by means of the Lagrange inversion: 


= ° + E 


r 


i! 


n=1 


f (x) r 


we obtain a formal series solving equation [TJ 

^ V 

x = a + > — 

' n 

n= 1 

Developing the derivative: 


7 \ Tl — 1 

a \ f e~ 


da; 


-6 


= a+V-r 

' n! 


‘n—1 

E 

,k=0 


n— 1 


n=l 

and remembering that: 

k / l 


i—l — k 


dx J \x — b 


(n + fc — 1)! / 1 


v da; / \x — b 
substituting, we found: 


°° l n S-l ( n - 1 


(n — 1)! \x — b 


n-\-k 


(-i)* 


= a +E^E 


n\ 

n= 1 0 


1 xfe(« + fc - 1 ) ! / 1 


l (-l) A 


(n — 1)! \a — b 


n-\-k 


ST' ^ ( n ~ !)• n — l — k na / i \ k i n ^ ~ 1)! ( 1 

a; = o+> — > —-rTTTT 71 fc e (—1) - ' 

' 7r! (n — 1 — fc)!/d 

n=l 0 v ' 


(n — 1)! \a — b 


n-\-k 


grouping terms not depending by k to the left of second summation sign: 


^ 1 fnle a \ nn ^ (n — 1 + k)\ ( -1 

= o+> — - r > T---- ' 

< nln. \ n — h In — 1 — 


n\n \a — bJ ^ (n — 1 — fc)!A;! \n(a — b) 
Remembering the definition of Bessel polynomials (see 0) : 


Bn(z)=J2 


(n + k)\ (z 


k—0 


(n — k)\k\ 


id) 


we can rewrite the solution in a more compact form: 
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( 2 ) 


E l / ILLC 

—i— ( - 

nln \a — 


nle a 


B n -1 


-2 


n{a — b) 


Symmetrically we can obtain another solution series by swap of parameters 
a and b: 


= *+£ 


n— 1 


nle° 


nln \ a 


-6 


B n -1 


+2 


i(a — b) 


( 3 ) 


2 A novel representation for Bessel polynomials 

It is interesting to observe that, in the development of above seen calculations, 
we have encountered a novel representation for Bessel polynomials: 


B n - 1 I —) =e~ x x n (-?- 
x \dx 


n— 1 r 


= E 

k—0 


(n — 1 + k)\ ( 1 

(n — 1 — fc)!(fc)! \ x 


( 4 ) 


The first cases for n=0..5, of the polynomials of typee X x n+1 (iiryy-e: 

d 


dx. 



e x " 

j 

X 1 


= 1 





x 


12 




x 


60 

x 2 


120 


x 


3 



20 180 840 1680 

x + X 2 x 3 X 4 


(d y 

e x 

= l-“ 

X 

420 

3360 

1 

15120 

30240 

\dx J 

_x 6 _ 

^ 2~ 

X z 

X 6 

X 4 

X 5 


whereas the classical Rodrigues representation for Bessel polynomials gives 
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( d 

6 * X (dx 


e*x 2 * 0 


= (-l)°l 


1 / d V r 

dx 


— 2*1 

e*x 


= (-l) 1 (l-2x) 




e* x 


2*2 


= (-1) 2 (l - 6x + 12x 2 ) 


_ _1 
e 33 




x 


3*2 = (_i)3 _ 12a; + 60a; 2 - 120x 3 ) 

= (-1) 4 (l - 20s + 180s 2 - 840s 3 + 1680s 4 ) 


_1 
e x 



(-1) 5 (1 - 30s + 420s 2 - 3360s 3 + 15120s 4 - 30240s 5 ) 


A direct derivation of the proposed representation of Bessel polynomials in 
4. obtained from the classical Rodrigues formula ([2]) will be provided in the 
next section. 


3 Rodrigues formula and transcendental equa¬ 
tions 


The appearance of orthogonal polynomials (i.e. Bessel polynomials) in solution 
of above seen transcendental equation is not surprising, indeed remembering the 
generalized Rodrigues formula my 


it can be placed in a form resembling the terms of the Lagrange inversion 
series, by suitable choose of terms Q and W: 


P n (x) 



fix) 


when f(x) is a product of suitable rational function and exponential in vari¬ 
able x, 

In following subsections, we will examine some examples of interpretation of 
Lagrange series terms as Rodrigues formulas. 
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3.1 Equation = le x and Laguerre polynomials 

We can remember p] for an employment of Laguerre polynomials for solution of 
a similar transcendental equation (another generalization of original Lambert 
W equation): 


= le x (5) 

x-t K ’ 

which can be formally solved by: 

* = <+£ («(* -»» ( 6 ) 

n—1 

where L (x) denotes the generalized Laguerre polynomials (inpQ, an equiva¬ 
lent formula based on first derivative of Laguerre polynomials was found). 

In the mentioned paper an non-trivial induction proof was given, but a 
more straightforward derivation can be obtained if we recognize that the above 
equation can be rewritten as: 


x — s = le x (x — t) (7) 

that can be solved by Lagrange inversion, obtaining for n-th term in the 
series (having placed: f(x ) = e x (x — t)): 

( j \ n—l 

-J e™{x-t) n (8) 

remembering the definition of Laguerre polynomials, the equation ([21) follows. 

3.2 Equation x = a + /e -:r2 / 2 and Hermite polynomials 

The equation: 

x = a + le~ x2/2 (9) 

leads to Lagrange term : 

Pn{x) = {iL) n le ~ x2/2 (10) 

remembering Rodrigues formula for Hermite polynomial (see [3]): 

H n (x) = e* 2 / 2 (J-'j e- 2 / 2 (11) 

we can rewrite the Lagrange series term (1101) by means of Rodrigues formula 
as: 


Pn{x) = e nx 2 H n _i (y/nx) 


( 12 ) 
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therefore the complete solution of equation ([!)]): 


^ in 

x = a + — e na / 2 H n -1 (y/na) (13) 

z —' n! 

n= 1 

3.3 Equation x = a + le eX and Touchard polynomials 

The equation: 


x = a + le e 

leads to Lagrange term: 



remembering Rodrigues formula for Touchard polynomials: 



therefore: 

P n (x) = e e T n -1 (x + login)) 
and the formal solution of (1 l ip : 

00 j^n 

x = a + 'S~ ' — e e “r„_i (a + log(n)) 

Z ' n! 


(14) 


(15) 


(16) 


3.4 Equation x = a + /x 2 e ^and Bessel polynomials 

The Bessel polynomials can be generated by the following Rodrigues formula 
(see [2]): 


B n {x) = 2~ n e% 



(17) 


An employment the Rodrigues formula for Bessel polynomials for solving! 
seems natural: 


, o 

x = a + lx e x 

therefore the formal n-th term in the Lagrange solution series is: 



1 Replacing x > 4 this equation becomes x = ax 2 + le 2x , having same form of [T] 


(18) 


(19) 
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that can be formally written in terms of Bessel polynomials as follows: f B n (.x)2 n e idx. 
Therefore the complete solution of equation H8l becomes: 


a; = a+y^ — / B n (a)2 n e “ dx 

, n' J 

n—1 

For the case in study equation 0, we have to place: 


( 20 ) 


f(x) = 


x — b 


3.5 Equation (x — a)(x — b) = le x and Bessel polynomials in 
novel Rodrigues-like form 

In previous section, starting from Lagrange series, we have found a novel Rodrigues- 
like representation of Bessel polynomials: 


-2 

x 


B n -i ( — ) = ( — 2)" x n e~ x 


Mr 1 

e x " 

\dx ) 

x n _ 


We now will provide an alternative proof; starting from the equation (6) we 
change x->2/x: 


B n {—)=e~* (x^ 
x J \ dx 


-)2n 


( 21 ) 

( 22 ) 


Now we will make use of the following identity (we will provide a proof in 
the next section): 


d_ 

d-z 


= x 


dx 


= x 


n+1 


„n—l 


so we can rewrite equation (ED: 


-2 

X 


B„ I — ) = e~ x x n+1 


r n+l 


(23) 


(24) 


that matches (0]), after rescaling x x/n. 


4 Proof of identity ( x 2 D) n = x n+l D n x n 1 by induc¬ 
tion 


Now we will provide a proof for identity : 


= {x 2 D) n = x n+1 D n x n ~ 1 

n cr. / 


(25) 
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where we have placed: D = -^ in order to save space. 

We start observing that, for n=l the following is trivial that: 

(x 2 D) = x n+1 D n x n ~\n = 1 (26) 

Therefore for n=l the identity holds. For induction we proof that, if the 
identity holds for n, then it holds also for n+1. Applying the operator (x 2 D^j 
to both members of (1251) : 

(x 2 D) ( x 2 D) n = ( x 2 D ) (x n+1 D n x n ~ 1 ) (27) 

(. x 2 D) n+1 = x 2 Dxx n D n x n - 1 = x 2 Dx n xD n x n ~ 1 (28) 

remembering the following commutation rules: 

Dx n . -x n D = nx™- 1 (29) 

D n x . —xD n = nD n ~ l (30) 

we can write: 

(x 2 D) (x 2 D) n = x 2 [Dx n .] [xD n ] x n ~ x = x 2 [x n D + nx "- 1 ] [D n x . -nD 71 - 1 ] x 71 " 1 

(31) 

= x n+2 D n+1 x n . +x 2 [nx n ~ 1 D n x. -x n DnD n ~ 1 - nx n - 1 nD n ~ 1 ] x n ~ 1 (32) 

Now we can show that the term in square parentheses is null, indeed applying 
above seen commutation rules: 


_ rr ,n-\-2 


= X 


X 


l .+X 2 [: 


nx^xD^ 


+ nx n ~ 1 nD n ~ 1 - nx n D n - nx^nD 11 - 1 ] 


(33) 


= x n+2 D n+1 x n .+x 2 [nx n D n - nx n D n + nx^nD* 1 - 1 - nx n - 1 nD n ~ 1 ] x n ~ x = x n+2 D n+1 x n . 

(34) 

Therefore: 


(x 2 D) n+1 = (x 2 D) (x 2 D) n = x n+2 D 


n+l x n ' 


(35) 


5 Concluding remarks 

Employing the Lagrange inverting series, a solution of the transcendental equa¬ 
tion (x — a)(x — b) = Ze x (that can be considered a quadratic generalization of 
the equation defining Lambert W-function) has been found in terms of Bessel 
orthogonal polynomials. 
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Once again (see ID), a transcendental equation can be formally solved by 
means of classic orthogonal polynomials, suggesting a link between Rodrigues 
formulas and the terms of Lagrange series. Moreover a novel representation 
for Bessel polynomials B n (^r) = e~ x x n+1 (^) [yf+rjis shown. Finally, an 

— -£r) = ( x 2 D ) = 

x n+1 D n x" _1 . for the reciprocal differentiation has been found. 

Further investigations on radius of convergence and numerical feasibility of 
above seen series (possibly by series acceleration techniques) are required. 

Moreover the occurrence of classical orthogonal polynomials in solution of 
mixed exponential/polynomial transcendental equations suggests to investigate 
for more general series of hypergeometric functions generated by Lagrange in¬ 
version of transcendental equations. 
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